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Abstract 

Recently we proposed that K fields, that is, fields with a non- 
standard kinetic term, may provide a mechanism for the generation 
of thick branes, based on the following observations. Firstly, K field 
theories allow for soliton solutions with compact support, i.e., com- 
pactons. Compactons in 1+1 dimensions may give rise to topological 
defects of the domain wall type and with finite thickness in higher 
dimensions. Secondly, propagation of linear perturbations is confined 
inside the compacton domain wall. Further, these linear perturba- 
tions inside the topological defect are of the standard type, in spite 
of the non-standard kinetic term. Thirdly, when gravity is taken into 
account, location of gravity in the sense of Randall-Sundrum works 
for these compacton domain walls provided that the backreaction of 
gravity does not destabilize the compacton domain wall. 
It is the purpose of the present paper to investigate in detail the exis- 
tence and stability of compacton domain walls in the full K field and 
gravity system, using both analytical and numerical methods. We find 
that the existence of the domain wall in the full system requires a cor- 
relation between the gravitational constant and the bulk cosmological 
constant. 
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1 Introduction 



Recently, the idea that the visible universe with its 3+1 dimensions is embed- 
ded into some higher dimensional space has received considerable attention. 
Consistency with observations requires that the propagation of matter and 
fields is restricted to this 3+1 dimensional subspace, at least for not too high 
energies. In this context, two slightly different scenarios have been studied. 
The subspace may either be of the topological defect type, in which case it 
has a finite, although probably very small, extension in the additional di- 
mensions. Some first proposals of this type have already been made in the 
eighties, see [TJ, [2J. In the last years, interest in this type of "topological 
defect" universes has increased significantly, and the name of "thick branes" 
has become customary for these objects. Some recent work may be found, 
e.g., in [3] - [7]. The other possibility is that the subspace is strictly 3+1 
dimensional, in which case it is known as a "three-brane" . Investigation of 
these three-branes started in the nineties, and the literature on this subject is 
too numerous to be quoted here. Some recent reviews, in which also further 
references may be found, are, e.g., [8] - [TO] . 

It is one of the crucial features of a "thick brane" cosmological model 
that a dynamical mechanism must exist which provides the confinement of 
all matter fields to the subspace. We recently proposed such a mechanism [TT] 
which uses as its main ingredients a scalar field theory with a non-standard 
kinetic term (K field theory) and the observation that topological defects 
with a compact support (compactons) exist in this theory [12]. For different 
aspects of K field theory, we refer, e.g., to [13] - [20], whereas the theory 
of compactons is discussed, e.g., in [21] - [23]. Interestingly, the complete 
suppression of the propagation of fields outside the support of the compacton 
is an automatic result of the model proposed in [TT]. Further, the propagation 
of linear perturbations inside the topological defect (i.e., inside the brane) is 
standard, in close similarity to the Kaluza-Klein reduction, in spite of the 
non-standard kinetic term. Specifically, there are no tachyons on the brane, 
and the evolution of linear perturbations is both unitary and causal. Inside 
the brane, the only remaining effect of the original K field theory resides in 
the values of the masses of the (Klein-Gordon type) linear fluctuation field. 

For the full theory, with gravitation coupled minimally to the scalar field, 
we showed in the same paper that, provided the compacton domain wall 
is not destabilized by the addition of gravity, bulk gravity solutions of the 
Randall-Sundrum type [21] (that is, localization of gravity on the brane) do 
exist. However, the analysis of the stability of the full infield plus gravity 
system and specifically of the existence of the compacton domain walls in 
the full system was not done in Ref. [TT]. It is the purpose of the present 
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paper to perform precisely this analysis. 

In Section 2, we briefly present the model together with the resulting field 
equations and their vacuum solution. In Section 3 we present some analytic 
investigations of the system without gravity for later convenience. In Section 
4 we turn to the analytic investigation of the full system with gravitation 
included. We perform a power series expansion about the boundary of the 
domain wall as well as about its center. Further, we establish some global 
properties of the resulting system of nonlinear evolution equations. The re- 
sult of this analysis is that the compacton domain wall does not exist for 
completely generic values of the two remaining parameters (the gravitational 
coupling k and the bulk cosmological constant A). Instead, the existence of 
the compacton requires a correlation between these two parameters, that is, 
for a given k the compacton exists if A takes a fixed k dependent value. In 
different words, the space of compacton solutions forms a one-dimensional 
curve in the k-A plane. In Section 5 we turn to a numerical invastigation 
of the full system. We determine the correlation between k and A (i.e., the 
above mentioned line in the k-A plane) to a high precision. Concretely, we 
fine-tune the value of k for a given A. Further we determine the compacton 
radius, that is, the extension of the domain wall in the transverse direction. 
In addition, we display figures providing some details of the numerical inves- 
tigation. In Section 6 we study the issue of stability of the compact domain 
walls established in the previous sections under linear fluctuations. Specifi- 
cally, we prove by a combination of analytical and numerical methods that 
the domain wall is stable under fluctuations of the scalar K field also in the 
full theory, with gravitational backreaction included. Section 7 contains a 
discussion of our results. 



where A is the cosmological constant , and X now includes the metric 



2 The model 



The action is (for details we refer to [TT] ) 





X = \g MN d M id N i. 



(2) 



Further, the potential is 




(3) 
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We will choose a 5D metric ansatz with a Minkowskian 4D slice, written in 
the form 

ds 2 = e~ A{y) (dt 2 - dx 2 ) - dy 2 . (4) 

We remind the reader that - up to coordinate transformations - this met- 
ric ansatz is the most general one having the Minkowski space symmetries 
ISO(3,l) on the brane, see e.g. [3|. For the scalar field we assume £ = 
Then the field equation for £ is 

- 8A y e y + = n\*ae - « 2 ) (5) 

and the independent components of the Einstein equations read 

\A yy -Al= l -[K + ^\\e-a 2 ) 2 ] (6) 

-Ayy = K £y. (7) 

At first sight one might think that the system is overdetermined because 
there are three ODEs for the two unknown functions £ and A, but the field 
equation for £ is not independent of the Einstein equations. Indeed, when 
one calculates the y derivative of Eq. (jSJ) and replaces A yy and A yyy with the 
help of Eq. (jZj) , one recovers Eq. (jSJ) • 

Remark: the derivation of the field equation for £ from the Einstein 
equations requires a division by k 2 and another division by £ y , therefore 
the field equation is a consequence of the Einstein equations provided that 
these two quantities are nonzero. In the vacuum sector, where £ y = 0, the 
field equation for £ is more restrictive and requires that £ = ±a, whereas the 
Einstein equations would allow for more general constant values of £. 

For later convenience, we rewrite the above system of equations in dimen- 
sionless quantities, by rescaling £ — > a£, V — > a 4 V, and z = Xy, A _1 A — > A, 
\k — > k. Further, we introduce the derivative C(z) = A z (z) as field variable, 
and get 

-\ct*+su = ae-i) (8) 

-C z -C 2 = l-A + K 2 (e-l) 2 (9) 



4 

\c z = « 2 e (io) 

This system still has the vacuum solution 
£ = ±1 = const. , A = ^-^\ z \ + const. =^ C — -^^sign(z) (11) 
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where 

A = -A>0 (12) 

is a positive constant. 

3 Non-gravitational case 

The system without gravitation has the field equation 

t% z = m 2 - 1) (is) 

with the compacton solution 

-1 ( Z - Zq )<-Z 

£(z) = { sinz -§< (z-zo) < f (14) 

Here Zq is the compacton center, which we shall frequently choose equal 
to z = 0. Further, the compacton reaches the vacuum values £ = ±1 at 
z± = Zq ± ^ . We will use the notation zq, z± for the center and boundary of 
the compacton, respectively, also in the case with gravitation. 

Although we know the explicit solution in the case without gravity, we 
want to study the behaviour of the solution near the lower boundary and the 
center for later convenience. We assume z = for the moment and introduce 
the new variable t = z — z^ = z + ^ (where £(i = 0) = C,(z = zJ) = —1 by 
assumption). The field equation remains unchanged, 

= - 1) (is) 

because the translation t = z — z_ is a symmetry. Now we insert the power 
series expansion 

f = -1 + b 2 t 2 + b 3 t 3 + ... (16) 
into the above equation (j!5p and find 



8bi, -2b 2 = (17) 



with the three solutions 



o.i.4 



42' 2 

where b 2 = corresponds to the vacuum solution and b 2 = \ corresponds 
to the compacton solution. b 2 = — | corresponds to a solution growing in 
absolute value for increasing t, which has infinite energy and is of no interest 
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in this context. Observe, however, that 6 2 = — \ together with £(t — 0) = 1 
would correspond to the anti-compacton. We remark that this degeneracy 
for the coefficient b 2 exists due to the nonlinearity in the kinetic term and 
would be absent in a theory with normal kinetic term. Once a choice for b 2 
has been made, the higher coefficients b n are determined uniquely by linear 
equations. For the choice b 2 — \, the higher coefficients just reproduce the 
power series of the function = — cost. 

Next, we study the power series expansion about the compacton center 
z = 0. By assumption, £(z — 0) = 0, therefore the expansion is 

£ = hz + b 2 z 2 + ... (19) 

If we assume that b\ 7^ 0, then it automatically follows that all even b n are 
zero, b n= 2m = 0. The odd b n are determined by linear equations depending 
on the arbitrary nonzero b\. For the choice b± — 1 we recover the power series 
expansion of the function £(z) = sinz. It is easy to see that finite energy 
requires, in fact, b\ = ±1. Indeed, a first integration of the field equation 
gives 

a = (e -if+k (20) 

where k is an integration constant. Finite energy requires k = which implies 
for £(zo) = that £ z (z ) = ±1. 

Finally, let us briefly recapitulate the necessary conditions for the exis- 
tence of the compacton. These conditions are that 

1) the field may be joined smoothly with smooth first derivative to its 
vacuum value. This holds because of the triple degeneracy of b 2 = 
0, |, — \ at the compacton boundaries z±. 

2) the field, which starts at the value = —1 at z_ and obeys the 
compacton boundary condition b 2 — |, evolves such that it takes the 
value £(zo) = at some point z . It may be checked easily that this 
property follows from the structure of the field equation. 

3) the field resulting from the evolution of point 2) is odd about the point 
z . This means that the evolution of £ for z > Zq will be the mirror 
image of the evolution for z < z and, consequently, will join the other 
vacuum, £ = 1, at z+ — 2 zq — Z— . This is automatically fulfilled in 
the case without gravitation, because a solution with initial condition 
£{zo) = 0, £ z ( z o) 7^ is always an odd function about z Q for arbitrary 
values of £ z (zo). 

We remark that points 1) and 2) in the above list continue to hold in the 
case with gravitation. On the other hand, point 3) will not hold in the generic 
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case. Instead, it requires a fine-tuning between the gravitational coupling k 
and the cosmological constant A. 



4 The case with gravity 

4.1 Expansion about the compacton boundary 

We again introduce the variable t = z — z_ where by assumption £(z = 
z_) = £(t = 0) = -1, £ t (t = 0) = 0, and C{t = 0) = -y^A/3. The Einstein 
equations are 

\c t -c* = -h+^te-if (21) 

\Ct = k 2 £, (22) 
and we insert the power series expansion 

e = -l + M 2 + M 3 + M 4 + ... (23) 



C = -yA/3 +cit + c 2 r + ... (24) 

For 6 2 we find again Equation ffTTj) with the three solutions fjlSj) . correspond- 
ing to vacuum, compacton and an infinite energy solution, respectively, like 
in the case without gravity. Choosing b<i = \ for the compacton, the higher 
coefficients are, again, determined uniquely by linear equations. Concretely 
we find 

1 FK , /A 1 /16A 1 \ . , 

&3 = -T5V 3 ' h = -y 3 420 (t5 3 + 2j (25) 

4 ~ "24 + 90 3 ' 6 ~ 2 • 3 4 ■ 5 3 ■ 7 V 3 J + 2 4 ■ 3 3 ■ 5 2 ■ 7 3 + 720' ( ' 

These coefficients now depend on the cosmological constant and go to their 
nongravitational values (the expansion coefficents of — cos t) in the limit A — > 
0, as they should. As 63 is negative, the cosmological constant tends to 
increase the size of the compacton. 

For the function C we find from Eq. (1221) that the coefficients C\ - C4 are 
identically zero. The two first nontrivial coefficients are 



Ak 2 8 
C5 "T5"' C6 "45 



k 2 JI/3. (27) 
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We remark that the first dependence of the 6j on the gravitational con- 
stant k 2 enters for the coefficient b 8 . The effect of the gravitational backreac- 
tion on the power series expansion at the compacton boundary is, therefore, 
negligible for not too large values of the gravitational constant. The con- 
tribution of k 2 to bs is positive, therefore the gravitational backreaction by 
itself tends to shrink the "compacton" radius. 

4.2 Approximate determination of the compacton ra- 
dius 

The fact that the expansion coefficients almost do not depend on the gravi- 
tational constant k 2 allows to determine the "compacton" radius \zq — z_| for 
a given cosmological constant A approximately. The idea is to use the power 
series expansion for £(t) to a certain order, and to determine the first zero of 
this polynomial for t > 0. We warn that for a fixed A and arbitrary k a true 
compacton does not exist. But the position of the first zero will almost not 
depend on k, therefore this position should determine the compacton radius 
to a good approximation whenever the compacton exists, at least for not too 
large values of A and k. Concretely, we will use the power series expansion 
up to sixth order, where k does not show up at all. 

Let us first check this proposal for the nongravitational case, where the 
exact compacton solution is known. This nongravitational solution is given 
in Eq. ( Tl4l) therefore the compacton radius p = \zq — z-\, variable t and 
function £(£) are 

P=f, t = l + z > £00 = -cost (28) 

Observe that all even derivatives of £(t) at the compacton center t — p are 
zero, because £ must be an odd function of z = t — p. Now we use the power 
series expansion of the function £(t) = — cos t about t = up to sixth order 
to determine the zero p approximately, see Eq. ( 1291) . 



m 


compacton radius p 


eit = z ) 


exact (— cost) 


tt/2 = 1.5708 





power series (0(t 6 )) 


1.5699 


-0.0208 



We find that the power series up to sixth order reproduces the exact value 
of the compacton radius to a precision of better than 10~ 3 . Further, we 
plot the value of the second derivative at the compacton center, £"(t = p), 
which must be zero for the exact solution. This value is off about 2% for the 
approximation. 
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cosmological constant A 



compacton radius p 





0.0001 
0.001 
0.0025 
0.005 

0.01 
0.025 

0.05 

0.1 



1.5699 
1.57145 
1.57481 
1.57767 
1.58092 
1.58554 
1.59488 
1.6058 
1.62233 



Table 1: Approximate determination of the compacton radius from the power 
series expansion 

Now we do the same calculation for the power series up to sixth order us- 
ing the coefficients 6j of the previous subsection for the case with gravitation 
and cosmological constant included, for different values of the cosmological 
constant. The result is presented in Table 1. It is clearly seen that the 
compacton radius increases with the cosmological constant. Later we shall 
compare these simple results with a numerical integration and will find that, 
indeed, they reproduce the compacton radius with a rather good precision. 
We repeat that the existence of the compacton requires a fine-tuning of the 
value of k, but that the value of the "compacton radius" calculated above is 
completely insensitive to the value of k. 

4.3 Expansion about the center 

Now we want to study a power series expansion about zq where by definition 
£(z = zq) = 0. Here we use the Einstein equations in the slightly rewritten 



where we choose the positive root for £ z , which describes the compacton (the 
negative root would give the anticompacton). As said, for a compacton £(z) 
must be an odd function about z . It may be checked easily that a necessary 



form 




(30) 




(31) 
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and sufficient condition for this is that C is equal to zero at z , too, C(z ) = 0. 
Therefore, the conditions at the center of a compacton are 

£(z ) = , C(z ) = 0. (32) 

These equations already demonstrate the possible troubles for the existence 
of a compacton. For suppose we start the evolution at the lower compacton 
boundary C,(z_) = — 1, C(z_) = — a/A/3 and for compacton boundary con- 
ditions (i.e. 6 2 = 1/2)- After some evolution £(z) will hit a point z such 
that £,{zq) = 0, but there is no reason that C(z) has evolved to the same 
value C(zo) = at z = zq. The condition C(zq) = is an additional condi- 
tion for the already completely determined initial value problem (evolution 
from the compacton boundary with compacton boundary conditions) and 
requires, therefore, a fine-tuning of the parameters k and A. 

A first conclusion which may be drawn easily is that A > for k ^ 
is a necessary condition for the existence of a compacton. For assume that 
we start the evolution of the fields at the center with conditions £(^o) = 0, 
C(z ) = for A = 0. The evolution equations are 

C z = \(C* + K\e-l?) (33) 

i z = -^(C 2 + ^(e-l) 2 ) 1 . (34) 



From these equations and the conditions C(z ) = 0, C z (z ) > it follows 
that C(z) is strictly positive for z > Zq, and from this it follows in turn 
that both C z and £ z are strictly positive for all z > Zq. Therefore, £ can 
never reach a point z + where £ z (z + ) = 0, which is the compacton boundary 
condition at the upper boundary of the compacton. The conclusion is that 
a compacton cannot exist for k ^ and A = 0. 

An upper bound can be easily found for A by inserting the power series 
expansion 

f = b lZ + b 3 z 3 + ... (35) 
C = c lZ + c 3 z 3 + ... (36) 

into the above evolution equations (we assume momentarily that z = 0). It 
immediately follows that 



W 2 A 
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which is real only provided that 

A < 3/t 2 . (38) 
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Some further qualitative conclusions may be drawn by investigating the evo- 
lution equations flHUj) . flHTj) . by assuming an evolution from the center z . 
The problem is that £ 2 must reach its value zero exactly at the same point 
z + where £ takes the value +1. When £ z has not yet reached zero at the point 
z + where £(z+) = +1, then it will never reach zero, because from this point 
onward both quantities C 2 and (£ 2 — l) 2 start to grow. This happens for too 
small A for a given k. On the other hand, if £ 2 takes its value zero before £ 
reaches the value one, then all the higher derivatives of £ at this point are 
singular, and the evolution equations break down at this point. This happens 
when A is too big for a given k. Therefore, for each given k^O, there exists 
precisely one fine tuned value of A such that £ z is zero exactly at the point 
z + where £ takes the value +1, and, consequently, the compacton exists only 
for this fine-tuned value. 

Remark: One may derive a single first order evolution equations for the 
"orbit" C(£), and this equation may be used to determine the correlation 
between A and k (but not for the compacton radius). Indeed, dividing (|30|) 
by (I3~T1) one easily derives the following equation, 



where 



Ct = /3 C 2 - 1 + a{i 2 - l) 2 4 (39) 



C=JjC, «- X , /?=3 — ) • (40) 



Between the compacton center and the upper compacton boundary (or, in- 
deed, between the lower and upper compacton boundary) both C z and £ 2 
are strictly monotonous, therefore this "orbit" equation indeed may describe 
the compacton. Integrating, e.g., from the center of the compacton to the 
upper boundary requires the "initial" condition (at the center) C(£ = 0) = 
and the condition (7(£ = 1) = 1. These two conditions for a first order ODE 
overdetermine the system and do not have a solution in general. Closer in- 
spection of the above equation ( !39l) again reveals that there exists exactly 
one f3 for given a (or one a for given (3) such that the compacton exists. 
Indeed, as C% always grows with growing a for fixed (3 (and always grows 
with growing ft for fixed a), independently of the values of C and £, the 
above statement follows. The numerical values of a and (3 (and, therefore, of 
A and k) such that the compacton exists may be easily determined numer- 
ically from Eq. ( |39l) and completely agree with the values determined from 
the numerical solution of the system of equations (1301) and ( |3TT) presented in 
Table 2 and Fig. 1. We present the result of a numerical integration for some 
fine-tuned values of k and A in the last figure (Fig. 15). Indeed, for these 
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fine-tuned values the function C(£) reaches the value one at £ = 1 to a high 
precision. 



5 Numerical calculations 



Here we present the numerical calculations for the Einstein equations ([H]), 
( fTUl) . First we display the results for the fine-tuned values of k for a given A, 
as well as the corresponding compacton radius, see Figures 1, 2 and Table 2. 
The fine-tuned value of k for a given A is defined by the condition that the 
numerical integration reproduces the compacton configuration. When k is 
too big (or, equivalently, A is too small), then £ x will never reach the value 
zero. When k is too small (or, equivalently, A is too big), then £ z reaches 
zero before £ reaches one. At this point a singularity forms and the numerical 
integration breaks down. The fine-tuned value for k is just the boundary case 
between these two cases. More precisely, the values for k 2 in Figure 1 and 
Table 2 are determined from a numerical integration which starts at the cen- 
ter of the "compacton" , that is, with initial values C = 0, £ = ("shooting 
from the center"). Within this numerical intergration, the displayed values 
for k, 2 are the smallest values for k 2 that can be found numerically such that 
a singularity does not form in the numerical integration. Strictly speaking, 
they are, therefore, upper bounds for the true fine-tuned values of k 2 , but 
they are determined to a high precision. 

In Table 2, we provide two values for the compacton radius. The first, p c , 
is determined by a shooting from the center, like k 2 itself. The advantage of 
this value is that it provides a strict lower bound on the compacton radius, 
because the "compacton" radius strictly decreases with increasing k 2 for a 
shooting from the center (here the "compacton" radius is defined as the point 
z + where £(-2+) = 1, which is at the same time the point where £ 2 has a local 
minimum), and the used value for k 2 is an upper bound. The disadvantage of 
p c is that its value is extremely sensitive to a precise determination of k 2 . An 
error in the seventh of eighth digit for k 2 , for instance, frequently translates 
into an error in the third digit for p c . Therefore, we determine the compacton 
radius in Table 2 also from a numerical integration which starts at the lower 
boundary of the compacton with compacton boundary conditions ( "shooting 
from the lower boundary"). We know already from the power series analysis 
of the last section that the so determined "compacton" radius pb almost does 
not depend on the exact value of n 2 (here the "compacton" radius is defined 
as the point t = z where £(£ = zq) = 0). Therefore, we expect the so 
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0.1 


0.2544852246 


1.59367 


1.60259 


0.2 


0.3732506647 


1.60912 


1.61565 


0.3 


0.4698973846 


1.62347 


1.62565 


0.5 


0.6332261243 


1.63673 


1.64145 


0.7 


0.7752995199 


1.65032 


1.65425 


1.0 


0.9666302322 


1.66660 


1.67030 


1.5 


1.252723578 


1.68682 


1.69224 


2.0 


1.514998322 


1.70648 


1.71058 


3.0 


1.999550771 


1.73613 


1.74101 


5.0 


2.886273449 


1.78229 


1.78828 



Table 2: Numerical determination of the fine-tuned value for k 2 and of the 
corresponding compacton radius. The displayed value of k 2 for a given A is 
the lowest value of k 2 that can be found numerically by a shooting from the 
center, such that a singularity does not form in the numerical integration. 
Therefore, the true value of k 2 could be slightly smaller, and the value in 
the table is, in fact, an upper bound. Correspondingly, the value of the 
compacton radius p c , which is also determined by a shooting from the center, 
is a lower bound. The values for k 2 are shown to a very high precision (with 
many digits), because the determination of p c depends very sensitively on 
a precise determination of k. Finally, we display the compacton radius p b 
which is determined by a shooting from the boundary and depends on k 2 
only very weakly. 
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determined value pb to better reproduce the true compacton radius than p c . 
We find, for instance, that p& obeys pb > (tt/2) for all A, which we know 
must be true at least for small A (see the power series expansion analysis of 
the last section). The disadvantage of pb is that it no longer provides a strict 
lower bound, therefore we display both values. 

Next, we present figures with some details from the numerical integration. 
We do the numerical integration in two different ways. First, we start at the 
center where both £ and C are zero ("shooting from the center"), see Figures 
3-8. We fix A = 0.5 and choose for k the fine-tuned value as well as one 
slightly smaller and one slightly bigger value. The results for other values of 
A are completely equivalent, and we choose the value of 0.5 for A just for 
convenience. 

Then, we present the corresponding figures which one obtains for the 
numerical integration which starts at the lower boundary of the compacton 
(more precisely, very near the lower boundary, where the compacton bound- 
ary condition b 2 = 1/2 is taken into account), and integrates all the way up to 
the upper boundary ("shooting from the lower boundary"), see Figures 9-14. 
We remark that from these figures one may see explicitly that the "com- 
pacton radius" pb determined from the shooting from the lower boundary 
depends only weakly on k 2 , and therefore, approximates the true compacton 
radius rather well even if k 2 is not exactly equal to its fine-tuned value. In- 
deed, from Figures 9, 11, and 13 we find the following values for p b , given 
A = 0.5 and the values for k 2 indicated in the table, 



A 


K 2 


Pb 


0.5 


0.633226 


1.64145 


0.5 


0.608400 


1.64236 


0.5 


0.640000 


1.64121 



(41) 



where the first value k 2 = 0.633226 is the fine-tuned value. We find that 
when k 2 has a relative deviation of about 10~ 2 from its fine-tuned value, 
then the relative deviation of pb is about 2 ■ 10~ 5 , i.e., quite small, indeed. 

Finally, we present the result of a numerical integration of the "orbit" 
equation ( 1391) for some fine-tuned values of k and A, with initial condition 
C(£ = 0) = 0- F° r the fine-tuned values C reaches the value C(£ = 1) = 1 
to a high precision. 

We conclude that the numerical analysis completely confirms the qual- 
itative discussion of the last section, both for an evolution from the center 
and an evolution from the boundary of the compacton. 
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6 Linear Stability 



In this section we want to discuss the issue of stability of the compacton 
domain wall under linear fluctuations. We shall not further comment on the 
issue of stability under gravitational fluctuations. Stability under fluctuations 
of the gravitational field for a whole class of thick brane models has been 
extensively discussed in Ref. [3J, and we demonstrated already in Ref. [TT] 
that our compact domain wall solutions belong to that class and, therefore, 
the analysis of Ref. [3j applies. 

On the other hand, we want to study linear stability under fluctuations 
of the scalar field for the system with gravitational backreaction (for the 
compacton without gravity, linear stability has already beed demonstrated 
in Ref. [H]). It turns out that stability of the full system with gravity 
may indeed be demonstrated by a combination of analytical and numerical 
methods, at least for not too large values of A (and, consequently, also for not 
too large values of k). The restriction to not too large values of A probably 
does not mean that stability does not hold for large values, but just indicates 
that an approximation used in our proof becomes too crude for large values 
of A. 

The linear fluctuation equation for the scalar field is 

#7„ + 2(U» ~ MDVz ~ (3£ 2 - = \e A e z O V (42) 

where □ is the four-dimensional Minkowski space wave operator, and A and 
£ are the compact brane solutions of Sections 4, 5. Further, 77 is the linear 
fluctuation. Outside the compacton region, where the fields £ and A take 
their vacuum values (specifically, £ z = 0), the only possible solution is 77 = 0, 
like in the non-gravitational case. For the case inside the compacton region, 
we use the separation of variable ansatz 77 — > rj(z)4>{x) and get 

D(f) = -oo 2 (j) (43) 

(the Klein-Gordon equation) and 

Vzz + 2(|^ - A z )r] z - ^-^V = -wVtj (44) 

where u 2 is the separation constant. After the transformation rj = j-fj, 
Equation (jUJ) can be brought into the Schroedinger equation form 

- Vzz + U(z)fj = u 2 e A f] (45) 
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where 

U{z) = A 2 Z -A ZZ + ^- 2A Z ^ + ( 10) 



where 



We now eliminate all higher derivatives terms (first and higher derivatives of 
£, second derivative of A) with the help of the field equations for £ and A. 
Further, we again use the notation C(z) = A z , C = and find 

9 27 3 V 3 s 
- 2^^£ + 4< 2 r-i - 2K 3 {1 ~ r e) \ 2 T- 1 i (47) 

r = ^(C 2 -l) + K 2 (£ 2 -l) 2 . (48) 

Stability requires u 2 > 0. This, in turn, follows from the positive semi- 
definiteness of the above Schroedinger operator 

H = -d 2 + U (49) 

because from (77 1 H \ rj) > and the obvious (rj\e A \rj) > it follows that uj 2 > 0. 

Next, let us study some properties of the potential U. U tends to plus 
infinity at the compacton boundaries z = ±p (here we assume that the 
compacton center is at z — 0). This follows from the fact that T goes to zero 
at the compacton boundaries, and from the fact that 

z^± P r 

Further, U takes its only minimum at the center z — 0, where £(z = 0) = 
and (7(2 = 0) = 0. The value at the minimum is 

^^ = °) = -50-(^ < 50 > 

The fact that U has an extremum at the center follows from the symmetry 
of U under the reflection z — > — whereas the fact that this is the only 
extremum (minimum) of U is not so easy to prove, because both the com- 
pacton solutions £ and C and the correlation between A and k are needed 
to demonstrate it. Here we restrict ourselves to a numerical demonstration. 
The potential U is plotted for the three values A = 10~ 3 , 0.5, 6.0 in Figure 16. 
The behaviour described above is clearly seen in these figures. It follows that 
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the Schroedinger operator with potential U has a discrete, non-degenerate 
spectrum of eigenvalues 

Hi> n = E n il> n , E 1 <E 2 <E 3 <... (51) 

and we want to demonstrate that E 1 = 0, E n > for n = 2, 3, — Observe 
that we have to expect one zero eigenvalue because of the Goldstone mode 
of the translational symmetry z — > z + a. 

The idea now is to find a new potential U which is less or equal to the 
potential U for all z in the compacton region, 

U(z)<U(z) , ze[-p,p] (52) 

and to study the spectral problem 

H$ n = (-& t + U)$ n = Ej> n , E 1 <E 2 <E 3 <... (53) 

for this new potential U. Stability is proved if the new potential has precisely 
one negative eigenvalue E\, whereas all the other eigenvalues E n ,n > 2 are 
positive. This follows from the inequality (which we prove below, see Eq. 

m) 

E n >E n V n (54) 

because then E 2 > E 2 > and, consequently, only the lowest eigenvalue 
Ei can be less or equal to zero. Finally, we know that it must be equal 
to zero, Ei = 0, because we know that a zero egenvalue must exist (the 
above-mentioned Goldstone mode of the translational symmetry z — > z + a) . 
Therefore, the Schroedinger operator H cannot have a negative eigenvalue, 
which is what we wanted to prove. 

Concretely, we choose the infinite square well potential for U, which is 
plus infinity at the compacton boundaries, and equal to Uq in the interior. 
The spectrum for this potential is well-known and is given by the eigenvalues 

where L is the width of the square well. In our case ^- = 1 and L = 2p and, 
therefore, 



7T 



2 



E n = —n 2 + U . (56) 
Stability requires E x < and E 2 > or 



n 2 



E 1 = — + U < (57) 

7T 2 

E 2 = — + U > (58) 
P 2 
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where Uq is given in Eq. floTJ]) and for p, A and k the values of Table 2 have 
to be used. 

Qualitatively we can observe at once from Table 2 that our stability proof 
cannot work for arbitrarily large A, because \Uq\ grows with growing A, 
whereas p~ 2 shrinks with growing A. The numerical results are shown in 
Figure 17 (for E{) and in Figure 18 (for E 2 ). We find that E\ < for all 
A as must hold because of the existence of the Goldstone mode. For E 2 we 
find that E 2 > at least for < A < 2, therefore we have proved stability 
for this parameter region. 

Remark: we could, in principle, determine the value of A where the transi- 
tion from positive to negative E 2 occurs with higher precision (it will happen 
very near to A = 3). But this would require the numerical solution of the 
compacton- gravity system for these values of A, which is quite demanding 
and does not produce important new insights, therefore we do not perform 
this calculation here. 

Remark: the fact that our stability proof does not work for high values 
of A does not necessarily mean that stability does not hold for these high 
values. It could just be the case that the approximation of the true potential 
U by the infinite square well U is worse for higher values of A. This is indeed 
indicated by the potentials U plotted in Figure 16. The higher the value of 
A, the steeper the valley in the center of the potential, and the worse the 
approximation by a flat line. 

Remark: the stability proof also works in the limit A — ► 0, k — > 0, that 
is, in the non-gravitational case. This may be seen most directly by using 
the non-gravitational compacton solution p — | , £ = sin z and the resulting 
potential U = —2 + 2 tan 2 z, which obviously may be bound by an infinite 
square wall potential with depth Uq = —2. Therefore, our stability proof 
also provides an alternative proof for the stability of the non-gravitational 
compacton of Ref . [12] , [TT] . 

Finally, we still have to prove the inequality (15^1) which we used in our 
stability proof. The proof only uses elementary quantum mechanics, but as 
we could not find it in any of the standard text books, we provide it here for 
completeness. 

Theorem: Given the two Schroedinger operators H = —d 2 + U{z) and 
H = -d 2 z + U(z) where U(z) > U V_z G (-p, p), U(z) > C/ V ze (-p, p), 
lim z -*± p U(z) = +oo and \\m z ^± p U(z) = +oo such that the spectrum of 
both H and H is purely discrete and non-degenerate, 

Hip n = E n ^p n , E 1 <E 2 <E Z < ... 

H?p n = E n $ n , E 1 <E 2 <E 3 < ... 
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and further it holds that 

U(z)>U(z) V ze(- P ,p) 

then the inequality 

E n > E n V n 

holds. 

Proof: We define the Schroedinger operator 

H{s) = -dl + U + s(U - U) (59) 
which has a purely discrete and non-degenerate spectrum 

#0)Vn(s) = E n (s)^ n (s) (60) 
for all s G [0, 1]. It follows that 

E n (s) = (Ms)\H(s)\Ms)) (61) 
and, after a derivative w.r.t. s, and using that 

((d/ds)ip n (s)\ip n (s)) + (il> n {s)\(d/ds)i/> n (s)) = 0, 

we further get 

^-E n {s) = {ip n {s)\{U - U)\i, n {s)) > 0, (62) 
from which the theorem follows immediately. QED. 

7 Discussion 

In this paper, we completed the investigation which we started in [TT| , study- 
ing in detail the existence of compact domain wall solutions in the full system 
of a quartic K field coupled to gravity. The result of this investigation is that 
such compact domain walls do exist, but their existence requires a correlation 
between the gravitational coupling k and the bulk cosmological constant A. 
Further, these compact domain walls are stable under linear perturbations. 
Firstly, we want to emphasize that the correlation between the gravitational 
coupling k and the bulk cosmological constant A we find in our system is qual- 
itatively different from a fine-tuning which is frequently found in three-brane 
models. There, a fine-tuning between the bulk cosmological constant and the 
brane tension (or some related parameters of the brane) is necessary in order 
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to achieve a vanishing effective four-dimensional cosmo logical constant. In 
our case, on the other hand, the full action density is always zero on-shell 
locally (i.e., for any local solution of the Einstein equations). The correlation 
is needed, instead, to guarantee the global existence of a solution with finite 
energy in the transverse direction in the non-gravitational sector alone, that 
is, the global existence of a topological defect or domain wall. This correla- 
tion is, in fact, a rather interesting observation by itself. One may see it as a 
drawback, requiring a "fine-tuning" of the cosmological constant for a given 
gravitational constant, or as an advantage, providing a "prediction" of the 
bulk cosmological constant from the value of the gravitational constant and 
the requirement of the existence of a thick brane universe. Quantitatively, 
the diameter of of the brane in the transverse direction grows with growing 
cosmological constant, and this growth may be seen as the net result of two 
competing influences. The gravitational constant by itself tends to shrink the 
brane size, whereas the cosmological constant tends to increase it. But these 
two constants are correlated, and the influence of the cosmological constant 
is stronger. The net effect is, therefore, an increase in the size of the brane. 

Another interesting point consists in the fact that we were able to go 
rather far in the analysis of the model, establishing both the existence and 
the stability of the compact domain wall solutions, as well as its numeri- 
cal properties to a high precision. Due to the inherent nonlinearity of the 
gravitational backreaction, such systems are usually quite difficult to anal- 
yse, and some simplifying assumptions (like replacing the "matter" fields 
by some effective energy momentum distribution) are frequently employed. 
In the present paper, on the other hand, we performed a full field theory 
calculation without any of these simplifying assumptions. 

Further, let us emphasize that the results of this paper should be easily 
generalizable to other models. In fact, all that is needed is the presence of a 
potential term for the scalar K field which allows for a vacuum degeneracy, 
and the presence of a generalized kinetic term which takes a certain non- 
standard form in the limit of low energy. Still, this latter condition is rather 
nontrivial, because the absence of the standard quadratic kinetic term at 
low energies is a necessary condition. In this sense these models are spe- 
cial, because the presence of the quadratic term is a typical feature of low 
energy effective theories. Under which conditions these special theories may 
be induced remains an open problem at the moment, which deserves further 
study. In any case, the fact that essential features of brane physics, like do- 
main walls of compact support, suppression of linear propagation outside the 
domain wall, and standard linear propagation inside the brane, are naturally 
provided by these models, makes it worth studying these special theories. 

Finally, we want to mention that recently infield theories have been used 
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within the context of brane physics for a slightly different purpose in [25]. 
In that paper it was shown that the distance between different (thin) branes 
may be stabilized if the bulk scalar field has a nonstandard kinetic term. 
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Figure 1: The fine-tuned values of k and A, for which a compacton exists, in 
the A-k 2 plane, for some selected values. 
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Figure 2: The compacton radius p for some selected values of A (and for the 
corresponding fine-tuned values of k 2 such that the compacton exists), in the 
A-p plane. 
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Figure 3: Shooting from the center: for A = 0.5 and for the fine-tuned value 
k 2 = 0.633226 the functions £ and £' = £ 2 are shown. The compacton field 
£ takes its vacuum value £ = 1 exactly at the point where £ z = (the com- 
pacton boundary). To the right of the compacton boundary the numerical 
integration selects the compacton configuration instead of the (correct) vac- 
uum configuration, because numerically £ 2 never is exactly zero (because in 
that case the numerical integration would break down). 
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Figure 4: Shooting from the center: for A = 0.5 and for the fine-tuned 
value k 2 = 0.633226 (the same values as in Figure 3) the functions C and 
C = C z are shown. The field C takes its vacuum value C = a/A/3 exactly 
at the point where £ z = C z = (the compacton boundary). The curve is 
very flat near the compacton boundary, because the first four derivatives 
of C at the boundary are zero. To the right of the compacton boundary 
the numerical integration selects the compacton configuration instead of the 
(correct) vacuum configuration, for the same reason as in Figure 3. 
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Figure 5: Shooting from the center: for A = 0.5 and for k 2 = 0.621335 
(slightly smaller than the fine-tuned value) the functions £ and £' = £ 2 are 
shown. The field £' reaches the value £' = before £ reaches its vacuum value 
£ = 1. Further, all higher derivatives of £ become singular at that point. 
Therefore, the numerical integration breaks down at this point. The numer- 
ical value of the point where the integration breaks down is z = 1.29653. 
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Figure 6: Shooting from the center: for A = 0.5 and k 2 = 0.621335 (slightly 
smaller than the fine-tuned value), like in Figure 4, the functions C and 
C = C z are shown. The field C reaches the value C = before £ and C 
reach their vacuum values. The numerical integration breaks down at this 
point because of the singularity in the integration for £, like in Figure 5. 
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Figure 7: Shooting from the center: for A = 0.5 and for k 2 = 0.633693 
(slightly bigger than the fine-tuned value) the functions £ and £' = £ 2 are 
shown. The field £' never reaches the value £' = 0, therefore the field £ never 
settles down at its vacuum value. 
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Figure 8: Shooting from the center: for A = 0.5 and for k 2 = 0.633693 
(slightly bigger than the fine-tuned value), like in Figure 7, the functions 
C and C = C z are shown. The field C never reaches the value C = 0, 
therefore the field C never settles down at its vacuum value (the fact that 
C does not reach zero is not obvious from the figure, because it takes very 
small values numerically near the region where £ ~ 1. But it may be checked 
easily by direct numerical calculation, or by amplifying that region in the 
figure) . 



29 




Figure 9: Shooting from the lower boundary: for A = 0.5 and for the fine- 
tuned value k 2 = 0.633226, like in Figure 3, the function £ is shown. The 
compacton field £ takes its opposite vacuum value £ = +1 exactly at the 
point where £ 2 = again (the upper compacton boundary). The graph of 
the function £ is exactly odd about the compacton center (where £ = 0). 
This compacton center t = Zq is at Zq = 1.6414529 



30 



c 



0,4- 



0,2- 







-0,2- 



-0,4-- 




Figure 10: Shooting from the lower boundary: for A = 0.5 and for the fine- 
tuned value k 2 = 0.633226 (the same values as in Figure 3) the function C 
is shown. The field C takes its opposite vacuum value C = +a/A/3 exactly 
at the point where £ z — C z = again (the upper compacton boundary). 
The curve is very flat near the compacton boundaries, because the first four 
derivatives of C at the boundaries are zero. The graph of the function C is 
exactly odd about the compacton center (where £ = C = 0). Numerically, 
this center is at t = 1.6414542, which agrees with the compacton center for 
£ in Figure 9 in the first six digits. 
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Figure 11: Shooting from the lower boundary: for A = 0.5 and for k 2 = 

0. 608400 (slightly smaller than the fine-tuned value), the function £ is shown. 
The field £' reaches the value £' = before £ reaches its vacuum value £ = 

1. Therefore, the numerical integration breaks down at this point. The 
numerical value of the point where the integration breaks down is t — 2.79688. 
The numerical value of the "compacton center" where = is t = 1.64236. 
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Figure 12: Shooting from the lower boundary: for A = 0.5 and k 2 = 0.608400 
(slightly smaller than the fine-tuned value), the function C is shown. The 
field C reaches the value C = before £ and C reach their vacuum values. 
The numerical integration breaks down at this point because of the singular- 
ity in the integration for £, like in Figure 11. The field C(t) takes the value 
zero at t = 1.67201, which is markedly different from the point where the 
field £ takes the value zero in Figure 11. 
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Figure 13: Shooting from the lower boundary: for A = 0.5 and for k 2 = 
0.640000 (slightly bigger than the fine-tuned value), the function £ is shown. 
The field £' never reaches the value £' = 0, therefore the field £ never settles 
down at its vacuum value. The point where £(t) = is at t — 1.641214. 
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Figure 14: Shooting from the lower boundary: for A = 0.5 and for k 2 = 
0.640000 (slightly bigger than the fine-tuned value), the function C is shown. 
The field C never reaches the value C — 0, therefore the field C never settles 
down at its vacuum value. The field C(t) takes the value zero at t = 1.63363, 
which is markedly different from the point where the field £ takes the value 
zero in Figure 13. 
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Figure 15: Integration of the orbit equation (1391) for some fine-tuned values 
of k and A. The values of the parameters and the resulting (7(£ = 1) are 
upper curve: A = 0.005 , k = 0.2303457790 , C{£ = 1) = 1.0000029, 
middle curve:_A = 1.5 , k= 1.119251347 , C(£ = 1) = 1.0000074, 
lower curve: A = 6.0 , k= 1.818434381 , , C((=l) = 0.9999953. 
So C reaches the compacton value C(£ = 1) = 1 to a high precision. 
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Figure 16: The potential U for the effective Schroedinger equation of Section 



6, for the values A = 0.001, 0.5, 6.0. 
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Figure 17: The first energy eigenvalue of the infinite square wall potential of 
Section 6. The eigenvalue is negative for all values of A 
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Figure 18: The second energy eigenvalue of the infinite square wall potential 
of Section 6. The eigenvalue is positive for sufficiently small values of A (at 
least till A = 2) 
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